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Abstract 

We compute 0(a s ) and 0(a™/3 { g _1 ) (BLM) corrections to the five structure functions relevant for 
b —> qlu decays and apply the results to the moments of a few distributions of phenomenological 
importance. We present compact analytic one-loop formulae for the structure functions, with 
proper subtraction of the soft divergence. 
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1 Introduction 


The experimental study of inclusive scmileptonic B decays has reached an unprecedented level 
of precision. The total branching ratio is now known to better than 2% [1], while various 
moments of the lepton energy and invariant hadronic mass distributions have been measured 
by the Delphi, BaBar, Belle, CLEO, and CDF Collaborations [2] with impressive accuracy. The 
first field of application of this wealth of data is the extraction of the CKM matrix element V c b, 
but the semileptonic distributions provide us with important information on the structure of 
the B meson and on the quark masses that can be employed in many applications. 

From the theoretical point of view, inclusive scmileptonic B decays can be described in the 
framework of the Operator Product Expansion (OPE) in inverse powers of the heavy b mass. 
The dominant term in this expansion represents the decay of a free bottom quark and the 
non-perturbative corrections to the parton model result are suppressed by at least two powers 
of the energy release, E re i « nib — m c [3]. The 0(1 /t?t|) [3, 4] and 0(l/mf) [5] corrections 
are known in terms of the B meson expectation values of a limited set of dimension five 
and six operators. These non-perturbative parameters can be extracted from the moments of 
the measured spectra of semileptonic and radiative inclusive decays, although complementary 
information is also provided by a number of heavy quark relations and sum rules (see e.g. 
[6]). Recent fits to the moments of the distributions lead to consistent results and allow for a 
determination of \V c b\ to better than 2% [1, 7, 8], with theoretical uncertainties starting to play 
a dominant role [9]. 

At the current level of experimental accuracy it becomes crucial to have perturbative cor¬ 
rections under control. In this respect, a lot of work has been done in the past 25 years. The 
one-loop corrections to a few differential distributions have been known for a long time [10, 11] 
and the correction to the total width is known in terms of polylogarithms since 1989 [12], An¬ 
alytic expressions for the 0(a s ) corrections to the moments of the hadron spectra have been 
computed in [13, 14], while numerical results for the leptonic and hadronic moments with a 
lower cut on the charged lepton energy can be found in [15] and [16, 14], However, the complete 
one-loop corrections to the five structure functions that enter the triple differential distribution 
for b —» civ have become available only recently [17, 18]. In the limit of massless leptons, just 
three of these structure functions are actually relevant, but their knowledge is necessary in 
order to compute the moments of a generic distribution with arbitrary experimental cuts and 
it is therefore phenomenologically important. At the level of physical observables, but not of 
structure functions, semileptonic b — > u decays can be obtained as the m c —> 0 limit of our 
results. In this case the expressions simplify considerably, and the complete structure functions 
are known since a few years [19]. 

Beyond the one-loop approximation, several results are already available. In particular, 
the O(al/3o) (i.e. the so-called BLM [20]) corrections to the total width have been computed 
in [21, 22], The equivalent corrections to the lepton energy spectrum can be found in [23], 
while the first results of a comprehensive analysis of hadronic moments have been reported 
in [18]. The BLM corrections are related to the running of the coupling constant in the loop 
contributions, and are practically most easily computed using one-loop diagrams calculated 
with a finite gluon mass [24], Starting from this result, it is moreover possible to compute and 
even resum higher order BLM contributions, O(a” +1 /3o) [22, 25]. Complete two-loop results 
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are known for specific kinematic points only (at zero, maximal, and intermediate recoil [26]), 
from which one can interpolate over the entire q 2 range and estimate the correction to the total 
width, which turns out to be well approximated by the BLM result. Recently, even three-loop 
O(a^) corrections have become available for a particular kinematic configuration, the extreme 
zero-recoil limit [27]. 

In this paper, we calculate analytically the complete one-loop corrections to the five structure 
functions relevant for the decay b —► civ, paying particular attention to the proper subtraction 
of the soft singularities, which is important for an accurate numerical evaluation. In this respect, 
our calculation differs from [17, 18]. In fact, our aim is to provide a reliable tool to compute 
the QCD corrections to the moments of B —> XJLv in the most general experimental setting, 
including the case i — r. Whenever possible, we give our results in analytic form; this in 
particular applies to the five structure functions and to a few moments with an arbitrary cut 
on the lepton energy. Our calculation represents an important check of previous results and 
provides new analytic formulae that allow for an easy and accurate numerical implementation. 
We find good agreement with previous results whenever a comparison is possible. 

Using the technique of refs. [24, 28], we also compute the O(a”/3o _1 ) corrections to the 
triple differential distribution and to its moments. To this end, we have performed analytically 
the calculation of the structure functions with a finite gluon mass, and integrate numerically 
the results over the gluon mass, applying the technique to the differential rate. We present 
numerical results for the one-loop and the BLM corrections in a few cases of practical interest. 

This paper is organized as follows: in Section 2 we set our notation, briefly describe the 
calculation of the next-to-leading order contributions, and discuss the subtraction of the soft 
singularity; we also report at the end of this section the complete expressions for the one-loop 
corrections to the structure functions. In Section 3 we describe the technique used to calculate 
the BLM corrections. Section 4 contains some reference numerical result and a comparison with 
the literature. Finally, we briefly summarize our results in Section 5. Details on the kinematics 
involved to next-to-leading order are given in Appendix A, while in Appendix B we illustrate 
the computation of the one-loop corrections. Appendix C gives analytic expressions for a few 
moments of the hadronic invariant mass and energy with an arbitrary cut on the lepton energy. 


2 The differential decay width to next-to-leading order 

The differential rate for the process 

b(p) -> £(pe) + v{Pu) + X c (r ) 


is given by 


dT 


G 2 F \V C 


cb\ 


Li*v(pe,Pf>) W^(p,q), 


dq 2 du dEi 16t^ 

where q = p( + pp, r = p — q, u = r 2 — m 2 , Eg is the charged lepton energy in the rest frame of 
the decaying b quark, and for massless neutrinos 


Pu) = Pe Pp ~ PePu 9^ +P u eP%- it^PtcJPvp- 


( 2 . 2 ) 
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Details on the kinematics involved are given in Appendix A. The decay rate is often expressed 
in terms of q 0 = (m b + q 2 — u — , m 2 c )/( 2 m b ), instead of u, which we have chosen for later 
convenience. The hadronic tensor W has the following, general decomposition, based on 
Lorentz covariance and CP symmetry: 

W^(p, q ) = -W 4 (q 2 , u) g pv + W 2 (q 2 , u) v»v v + iW 3 (q 2 , u)e^ pa v p q a 

+W 4 (q 2 , u) <f + W 5 (q 2 , u) (v p q u + v v q% (2.3) 


where v = p/m b , q = q/m b , and u = u/m b (in the following, we denote with a hat dimensionless 
quantities, normalized to the b quark mass). In the limit of massless leptons the terms propor¬ 
tional to IT 4 and W 5 do not contribute to the decay rate, since in that case q^L pu = q v L pv = 0. 
One finds 


m 2 b 


Q 2 W\(q 2 , u) 


2 (vp e ) 2 - 2 vptvq + y 


W 2 (q 2 ,u) + q 2 ( 2 vpi 


vq) W 3 (q 2 ,u). (2.4) 


The squared matrix element depends on the charged lepton energy E( only through vpi = E(. 
The leading contribution to the hadronic tensor W pi/ describes the process 


b(p) -> e(p r ) + v(p„) + e(p') 


(2.5) 


at the tree level. One finds 

Wftip, q) = [-Tj g pu + T 2 v p v v + %T 3 e^v a q p + T 4 q p q v + T 5 (q p v u + q u v p )\ tt S(u), (2.6) 

where 


T 1 (q 2 ) = -tu T 2 (q 2 ) — 4 T 3 (q 2 ) = 2 T 4 (q 2 ) = 0 T 5 (q 2 ) = -2, (2.7) 


and we have defined 


u = q 


P\ 


m : 
P =~2 


m 


2 • 


( 2 . 8 ) 


The hadronic tensor receives contributions of order a s from one-gluon emission at the tree 
level, 

b(p) -*• e(p t ) + D(p v ) + c(p') + g(k), (2.9) 


integrated over the charm-gluon phase space (in this case, r — p' + k), and from one-loop virtual 
corrections to the process in eq. (2.5): 


» r ,y<P. 9 ) = Wfi R (p,q) + Wfiyfaq). (2,10) 

The computation can be performed either by computing the relevant squared matrix elements 
and integrating over the corresponding phase space measures, or by taking the imaginary part 
of the bW* —>■ bW* forward amplitude. We calculated the structure functions in both ways; in 
Appendix B we give some details about the first method only. We specifically consider here only 
the weak decay structure functions and assume the corresponding kinematics, i.e. ml>q 2 > 0 
and q 0 > \fq 1 . 
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Let us begin with the contribution from real emission. We have performed the computation 
with the help of the algebraic manipulation programs Maxima [29] and Mathematica [30] for 
an arbitrary gluon mass A. This is necessary for two reasons: first, the gluon mass regulates the 
logarithmic infrared divergence that arises in the phase space integration from the soft part of 
the gluon spectrum, and second, the full dependence of the one-loop corrections on the gluon 
mass allows us to compute the order-(a:”/3o _1 ) corrections. The calculation is straightforward, 
and is already discussed in the literature in different contexts [11]; the only point that deserves 
a careful discussion is the isolation of the infrared-singular term, that is canceled by a similar 
contribution from virtual diagrams. Both for conceptual reasons and for the numerical imple¬ 
mentation, it is important that the final result is manifestly independent of the gluon mass in 
the nig —> 0 limit. 

While tree-level and virtual contributions are characterized by a S(u) that defines their 
kinematic structure (production of an on-shell charm quark without gluons), the contributions 
from real emission diagrams are defined for values of the variable u between 1 

U- = 2Xy/p + A 2 and u + — (1 — \J~<p) 2 — p. (2.11) 

In this range, the real emission contributions to the structure functions are regular as A —»• 0, 
but practical applications require an integration over u. In the physical limit A —^ 0, when 
the integration is performed close to the kinematic boundary u = U- = 0 the real emission 
contributions develop a log A singularity. As mentioned above, this singularity is canceled in 
physical quantities by a corresponding term originating from the virtual corrections. Since the 
singularity is localized at u — 0, it is possible to extract it and write it in the form of a term 
proportional to 5 (it) log A, making the cancellation of the infrared divergences explicit already 
before integration. 

This can be done as follows. The full result can be written as 

w Tm(P’ 9) = w m ,,<«(?, 9) + 9) (2.12) 

where we have isolated in W(^ Rsing (p, q) all the terms that give rise to a soft singularity when 
integrated over u. There are obviously different ways to perform the splitting between the two 
above contributions, that differ in IL/]/ R sing by terms that do not give rise to singularities upon 
integration over u. Let us first concentrate on W^ Rsing , which we define as 

W[.? )R , iJm (p,<l) = C F a, [i» 9 '“' + 4t,V + 2i€'“'"' , K a ®-2(K'T + Vr)] D(f, fi.A 2 ), (2.13) 
where Cp = 4/3, and 



ICj t + 1 f_|_ — t_ 
u y/Xf g 4 y/\b i + t_ 


r x 2 . (h + A 2 )(cb - u) ± V%y/X~ c 

£+ — A H-—--- 

2(w + p) 

X b = (Cj - u) 2 - 4 (u + p) 

X c = {u- A 2 ) 2 - 4A 2 p. 


P_ 4 ~ t- 
u 2 \fX b 


(2.14) 

(2.15) 

(2.16) 
(2.17) 


1 On general grounds, the structure functions are actually proportional to 9(u — it-). In the following, we 
omit all 9 functions from the explicit expressions, as it should not cause confusion. 
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In the limit A —► 0 


D(q 2 ,u, 0) = 
D 0 (q 2 ,u) = 


D 0 (q 2 ,u) 


u 


2 _ Cj Cj - U + y/\ u + 2p 

y/Xj, Co — U — y/\b U 4- p 


(2.18) 

(2.19) 


so that the integral over u is logarithmically divergent. The two above equations are sufficient 
to isolate the divergent term, but not to evaluate the finite contribution. The problem is that 
D(q 2 ,u, X) is not analytic when u = U- and A = 0: indeed, 


D(q 2 ,ii, A 2 ) = D s (q 2 ,ii, A 2 ) + 0(u - tt_, A) 

1 a i , VX, V» 2 -» 2 - 

D s (q 2 ,u, A 2 ) = — log- 7 =— , = -rg —2 

u ^ \/u 2 -u 2 _ -2 _l Vffi 


u 2 — u 2 _ 


u 2 — u 2 _ 


u 2 + 




( 2 . 20 ) 

( 2 . 21 ) 


4p 


where A° = d> 2 —4p. The integration over u becomes non-trivial when it starts from the endpoint 
u = U-. We now add and subtract to D(q 2 ,u, A 2 ) a term proportional to 5(u): 


D(q 2 ,u, A 2 ) = D(q 2 ,u,X 2 )-A(q 2 ,X)5(u) +A(q 2 ,X)S(u) 
where we have defined 


r u + 


A{q i A) = / duD s (q 2 ,u, A 2 ). 


( 2 . 22 ) 


(2.23) 


The u integration of the term in square brackets in eq. (2.22) is regular for A —» 0 over any 
range, while the soft singularity is entirely contained in A(q 2 , X). For instance, for a generic 
test function F(u) regular in u = 0, we have 


c u + ~ r u + 

du F(u)D(q 2 ,u, X 2 ) = / du 


F(u)D(q 2 , u, A 2 ) - F(0)D s (q 2 , u, A 2 ) + F(0) A(q 2 , A), 


(2.24) 

The integration range in eq. (2.24) is irrelevant and may be restricted by the shape of the test 
function. We have chosen it to coincide with the allowed kinematic range < u < u + so that 
for F(u) = 1 the integral in eq. (2.24) is the one which is relevant for the computation of the 
total rate. 

One can make notation more compact by defining 


D(q 2 , u, A 2 ) = [D(q 2 , u, A 2 )] + + A(q 2 , A) 5(u), 
where the plus prescription is 2 


ru+ 


du 


D(q 2 ,u,X 2 ) t F(u) = 


r u + 


du 


D(q 2 ,u , X 2 )F(u) - D s (q 2 , u, A 2 )F(0) 


(2.25) 


(2.26) 


2 This is a generalization of the standard ’phis’ distribution. 
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In the limit A —> 0, we have 


D(q ,u, A ) , = D 0 (q , u) - 

u 


where the distribution (1 /u)+ is dehned by 


J du (i') f(u) = f + ^[/(h)-/(0)]. 

Jo \uj+ Jo u 


(2.27) 


(2.28) 


As anticipated, the last term in eq. (2.25) has the same kinematic structure as the tree-level 
and virtual contributions. The integral A(q 2 , A) can be computed analytically for generic values 
of A; in the small- A limit we have 


A(q 2 , A) = ( 2 - - log log + J- 
' a 1 — a) 


u + 

A log (■—— + 1 + e(X 2 ), 

la l — o 


2 a 


Lio 


2 a 

a — 1 


Li? 


2 a 

a T 1 


(2.29) 


where a = \Jx\/C j. Finally we have 


Wising (p,q) = C F a s [iog^ + 4v^ + 2ie^v a q p -2(iAq" + v u qn\ [D(q 2 ,u, A 2 ) 
Cp Ot-s 


+ 


7T 


Wft(p,q)A(?,\). 


(2.30) 


We now consider the regular part of the real emission contribution. It can be written in 
terms of structure functions as in eq. (2.3): 


W tm,reg(Pi q) = C F a s {- Ri(q 2 , u) g^ u + R 2 (q 2 , u) v' l v v + iR 3 (q 2 , u) e^ a(5 v c 

+ R±(q 2 ,u) <f + R§(q 2 1 u) (v M <f + gV')}. 

The analytic expressions of the structure functions R, r for A = 0 are 
Ri(q 2 ,u) 


% 


(2.31) 


u (3m + 2p) + to (3u + 4p) | u 2 {oj — u) | 6 u 


4 (it + p ) 

4h 2 — (6u) + A b)u + 2X b A 

2\b\A\b 


2 (u + p) X b X b 
log t 


(2.32) 


R'2(q 2 ,u) = 


2 [am(8 — 9 u) — 8a) 3 + 15a) 2 h + 2u{u 2 — 8u + 8)] 64 (1 + a) — u) p 


A?, 


A;, 


+ 


am 


(a) — h) 3 2 (2a) 4 + a) 2 -u(7 + 4u) — 5d) 3 h — 9dm 2 + 2h 3 — dm 3 ) 


Aft (u + p) 

2 


A;, 2 (h + p) 


2A& + Aft(26h — 2u 2 + a)(8 + 3 it) + 16p) 


Xl\fXb 

—6 u (3u 2 + 2u(p — 5) — 12p — d)(3 + 4h + 3 p)J 


(2.33) 

(2.34) 


log r 
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RsiflU) = 


(Co — 2 )u — Co 2 — 4(a) + p) + Xb 


AbV^Ab 


logr 


P 


2(u + p)" 


3a) + p 


2 (a) + p) (u + p) 


8a) + 5a) 2 — 3a)w + 4p + 4d)p — 2-up 
Af, (a) + p) 


R4g 2 ,u) = 


8 (2a) 2 — 2u — 5am + 3 m 2 ) 64p u 2 (Co — u) 3 \3u(2u — 3a)) 

+ 2" + 


A h 


Xb Xb (u + p) 


X h 


(2.35) 

(2.36) 


+w — 2a) 


4 log r 2m (m 3 — 3a) 2 (l — w) — a) 3 + Co(u — 3 u 2 ) + 2u 2 ) 


R 5 (q 2 ,ii) = 


AfeV^ At 2 (m + p) 

2 (4a) 2 + 2a) 3 — 12 m — 2am(3 — 2u) — 5Co 2 u + 2 u 2 — it 3 ) 8 (8 + 2a) — m) p 


A/, 


A h 


u(Co — u) 3 (Co + u) 2a) 4 + a) 2 w(10 + u) — 4a) 3 w — 2um 2 (5 — w) — w 4 


2X b 2 (u + p) Ah (ft + p) 

Ah(16w— m 2 + a)(12+-u) + 8p) + 6u(8u — u 2 + a)(6+w) + 12p) 


(2.37) 


log^ 

A 2 \/Ah 


J 


where 


oo — u 


r = - -7-7y=- (2.38) 

a i — u — VA b 

We now turn to the virtual contribution to the differential width. We calculate the one-loop 
corrections to the decay of an on-shell b quark into an on-shell c quark and a virtual W boson, 
parameterizing the result in terms of pole masses. Some detail of the calculation is given in 
Appendix B. The result has the form 


Wfi y (p,q) = -C F a s {-V 1 (q 2 ,X)g^ + V 2 (q 2 ,X)v^ + xV 3 (q 2 ,X)e^v a q p 

C p OL 


+ w, a) rr + w, a) k<t + w - 


71 


(2.39) 
W^(p,q)V„(q 2 ,X), 


where 


ko(4 A) = - (1 + J 0 — 2 J — 2/x) 

Ah(<? 2 ,A) = —CoK — (Co 2 — Co — 2p)I x — (Co 2 — Cop — 2p)i y 
V2(<r, A) = —4 (a) + p + 1)4 — Co(I x + 4 
V3(g 2 , A) = 2A + 2(a) — l)-4 + 2(a) — p)4 
4 4 (<f,A) = 4(4-4) 

AA(4; A) = 2 I"(a) + 2)4 — (ca + 1)4 — (a) — p)4 


(2.40) 


and the integrals 4, J, etc. are defined in Appendix B, eqs. (B.8) and (B.13). In the limit 
A —> 0 we have 

i 1 , A 2 y 1 f 1 z + 1 “ P , 

J = -l-log—; A = 4 = 4^ log-^logp 

4P 2 2<4 2 q 2 


(2.41) 




where the dots stand for terms regular in the limit A —> 0. It cancels exactly the singularity we 
have found in eq. (2.29). 

We conclude this section by summarizing our results for the structure functions to order a s 
for massless gluon. We have found: 


wm\u) = 

Trj Ti{q 2 )5{u) 

a 

+ - 

S C F 

71 

W 2 ) 

S(q 2 ) ~ VM 2 \ 

)] 




a 

+ - 

sC F 

m 2 ) 

D 0 {q 2 ,u) 

) +Ri(q 2 ,u) 

}. (2.47) 




71 

\uj 

+ 

where we have defined 



A{q ' 






S(q 2 ) = 

lim 

Wo 

2 ,A)- 

Vo{q 2 A) • 


(2.48) 

The T,j are given in eq. 

(2.7), A in eq. 

(2.29), V, 

; in eq. 

. (2.40), Ri in 

eqs. (2.32-2.38) 

, and D 0 in 


eq. (2.19). 

Ref. [17] also reports the results of the calculation of the structure functions W t . However, 
a direct comparison of our results for IT 7 ,; with ref. [17] is possible only before subtraction of 
the soft singularity. We agree with all the real emission contributions of ref. [17], but the 
virtual contributions to W 3 , W 4 , W§ (in the notation of ref. [17]) are different from our results. 
The discrepancy, however, is such that our Vi,V 2 , V :i are reproduced correctly. Therefore, a 
difference should be observed when computing the decay rate into massive leptons only. 


3 Generalized BLM — 0( ol ™(3 q j ) corrections 

The full perturbative expansion for the total rate can be written as 


r 


r 0 


1 + ^( 0 ) — + 

71 



(3.1) 
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where 


r ° = ^p(l - 8p + 8p 3 - p 4 - 12p 2 lnp) (3.2) 

is the tree-level semilcptonic width, and the argument in the order-aq correction cq is the gluon 
mass, which is zero in the physical limit. The coefficients a n+ 2 are degree-(n + 1) polynomials 
in the number of light fermions rif , or equivalently in the first coefficient of the (3 function for 
a s , /3o = 11 - 2n//3: 

a n +2 = a^+ 2 1 Po +1 + O(Pq). (3.3) 

The dominant term in the large-A limit, usually called the (extended) BLM approxima¬ 
tion [20], can be computed in terms of the first perturbative coefficient ai(A 2 ) [24]; in the MS 
scheme for the coupling constant the result reads [9, 24] 


a n +2 ~ 


■A nk 


2 \k ( Tl + 1 

J [ 2k+ 1 


ro o /7\^ //^ 

_ iog n - 2fc — 

d A 2 S A 2 


ai(A 2 ) — «i(0)- 


yU 2 = m 2 e 5//3 . 


Jo A 2 A 2 [ ' LK 'n 2 + \ 2 \' ^ • 

We will now obtain an expression similar to eq. (3.4) for the BLM correction to the differential 
rate. As we have seen in the previous section, the order-aq correction to the total rate is obtained 
by summing the virtual-soft contribution proportional to 6 (u) and the real contribution: 

«i(A 2 ) = a7 s (A 2 ) + af(A 2 ). (3.6) 

After integration over the lepton energy Eg, with or without kinematic cuts, we have 

af (A 2 ) = I'* dq 2 (3.7) 

Jo dq 2 

af(A 2 ) = e(S - A) f tS ~ X)1 df r du da A2) (3.8) 

Jo J u— dq 2 du 


where we have defined <5 = mb ~ nv The integration bounds on u = r 2 — m 2 are 

■u_ = (m c + A) 2 — m 2 ; u + = (mb — xfq 2 ) 2 — m 2 . 


u- — (m c + A) 2 — m 2 ; u + — (mb — \jq 2 ) 2 — m 2 . (3.9) 

The u integration in the virtual-soft term is trivial, since it is proportional to <5(w). Also note 
that the real emission term vanishes when A > 5; this is implicitly taken into account in the 
definition of a j, which contains the soft contribution of the real emission term. We have 


j[”^log“- 2 ‘^Lf(A 2 ) + af(A 2 ) 


A 2 6 
r* 2 d\ 2 

Jo 'A 2 ' 
d\ 2 , 


A 2 

«] 

h 2 

r ‘ 2 

A 2 

Jo 

2fc 

h 2 


A 2 . 

-2/ 


+ a?^ 2 ) - n? + a ?(o)) 

da\ s (q 2 , A 2 ) da\ s (q 2 , 0) fi 2 
dq 2 dq 2 A 2 + fi 2 


2 f u + j daf(q 2 ,u, A 2 ) 

tq L dU d q Mu 

r u + A daf(q 2 ,u, 0) ji 2 

3 “ dfdki A 2 + /x 2 


(3.10) 
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Now we observe that 


rS 2 


dX 2 


with 

Thus, 

A-nk 


f°° d\ 2 r 52 


X 2 


dq 2 


r<5 2 


dq 2 


ru+ 


5—A) 2 

dq 2 

J 

f u + f 5 „ 

/ du= d,q~ 

'u - Jo 

ru+ 

/ du 
Jo J 

,A 2 

/ dX 2 , 

'0 

A 2 

= (^Vr 2 — m^j . 



da x L s (q 2 , A 2 ) da\ s (q 2 , 0) 

d 2 

log"- 2 * L 

dq 2 

dq 2 

X 2 + /i 2 _ 

r A2 dX 2 

daf(q 2 , xx, A 2 ) 

daf (q 2 , xx, 0) /x 2 

Jo X 2 

dq 2 du 

dq 2 du A 2 + ji 2 


(3.11) 

(3.12) 


rS 2 


dq 2 


du 


daf(q 2 ,u, 0) r°° dX 2 /x 2 n _ 2k fi 2 
3 dg 2 dxx A 2 A 2 A 2 + fi 2 ° S A 2 ‘ 


log' 


i—2fc / i 


A 2 

(3.13) 


The A 2 integral in the last term can be expressed in terms of polylogarithms; for example, 


r°° dX 2 /i 2 

/a 2 A 2 A 2 + /i 2 


= log 


A 2 + jj 2 


A 2 


r°° dX 2 g 2 


d 


/a 2 A 2 A 2 + /x 2 & A 2 


log 


= Li, - 


d_ 

A 2 


- log — log 
d~ 


/x 2 + A 2 


A 2 


(3.14) 


From eq. (3.13), we can read off an expression for the BLM corrections to the differential rate: 


dA-nk 

dq 2 du 


f°° dX 2 


A 2 


h(xx) 


da{° (g , A 2 ) dcq b (xy,0) ji 


dq 2 


dq 2 


X 2 +/I 2 


log 


n—2k d 

A 2 


r A2 dX 2 

dalf(q 2 , xx, A 2 ) 

daf (q 2 , xx, 0) 

d 2 

Jo A 2 

dq 2 du 

dq 2 du 

A 2 + /i 2 _ 


log"- 2 * ^ 


daf(q 2 ,u, 0) [°° dX 2 /i 2 


dq 2 du J a 2 A 2 A 2 + /x 2 


log' 


i—2 k d~ 

A 2 ' 


(3.15) 


Eq. (3.15) can be used for a numerical computation of the BLM coefficients, since each integral 
is separately convergent. Any kind of kinematic cuts can be implemented using appropriate 
step-functions in the integrands. 

All the one-loop results of Section 2, as well as the above expressions, have been derived in 
the on-shell scheme for the b and c quark masses. The use of pole masses is known to induce 
large higher order perturbative corrections and is not adequate to precision studies. In most 
practical B physics applications (see e.g. [31]) it is therefore appropriate to employ well-defined 
short distance mass parameters, such as the so-called kinetic masses, namely heavy quark 
parameters which are renormalized a la Wilson and depend explicitly on a ‘hard’ normalization 
scale g which is chosen close to 1 GeV [25, 32], The change of scheme required to express 
0 (ck” +1 /9q) corrections in terms of short-distance parameters is elementary, once the relation 
between the latter and their on-shell counterpart is known at the very same order. 
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n 

LW 

n 

n 

n 

LIT 

n 

IT 

l. 

-1.777588 

-1.917079 

-2.827995 

1 

0.307202 

-0.551243 

-0.617666 

-0.938365 

2 

0.103000 

-0.187687 

-0.217405 

-0.338429 

3 

0.036524 

-0.067828 

-0.080868 

-0.128695 


Tabic 1: Leptonic moments without Eg cut. 


4 Numerical results 


The primary held of application of the results described in the previous sections is the calcu¬ 
lation of the perturbative contributions to the first few moments of various differential decay 
distributions. The moments most frequently employed are those of the charged lepton energy 
spectrum, and of the invariant mass of the hadronic system in the final state [2]. As experimen¬ 
tal cuts are usually applied, the moments of phenomenological interest are actually truncated. 
A typical case is a lower cut on the charged lepton energy. At the one-loop level, a fully analytic 
approach is possible even with a cut on the lepton energy, but the integration over the phase 
space is a tedious task. We have performed the phase space integration analytically for several 
of the partonic moments that represent the building blocks of the first three hadronic mass mo¬ 
ments (see [13, 31]), in the case of an arbitrary cut on the lepton energy. A few representative 
examples of the results are reproduced in Appendix C. However, in general, a strategy for an 
accurate and fast numerical evaluation is necessary, especially when computing higher order 
corrections in the BLM approximation. 

We have implemented these calculations in a FORTRAN code, which is available from the 
authors upon request. The subtraction of soft singularities described in Section 2 allows us to 
compute a generic moment and to apply arbitrary kinematic cuts. 

In this section we provide a few tables of reference numbers obtained with our code for 
massless leptons; they have been checked in a number of ways, both analytical and numerical. 
All results are in the pole mass scheme. 

In particular, we focus here on the leptonic moments, 


L n = — / dq 2 du dEp E 


d 3 T 


T 0 ■’ 

= 4 0) + 


dq 2 du dEp 


a sPo T ( 2 ) , 


I T (i) l T (2) , r (3) , 

i o J >1 i o 


n 


7T Z 


7T 


(4.1) 


and on the hadronic moments of the form 


Hij = F 
J- o 


dq 2 du dEp u 1 Eq 


d 3 T 


dq 2 du dEp 

,3/32 


/A°) , rr(l) , a sPo tt( 2) a s Po rr(3) 

11 ij ' 11 ij ' o 11 ij ' q 11 ij ' •** 

•y /jy J yjY ^ J J 


(4,2) 


where T 0 is the total tree-level width of eq. (3.2), and E 0 = (ml + r 2 — q 2 )/2mb is the energy 
of the hadronic system at the partonic level in the b quark rest frame. We employ the on-shell 
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n 

LW 

n 

LW 

n 

LW 

n 

Id 3 ) 

n 

IT 

0.814810 

-1.439726 

-1.600023 

-2.407834 

1 

0.277642 

-0.497511 

-0.567236 

-0.870422 

2 

0.097933 

-0.178501 

-0.208715 

-0.326798 

3 

0.035615 

-0.066185 

-0.079299 

-0.126665 


Table 2: Leptonic moments with Eg > 1 GeV. 


i 

j 

H m 

H [1) 

*7 

h (2) 

i. i 

h (3) 

IT 

0 

1. 

-1.777615 

-1.917079 

-2.827995 

0 

1 

0.422009 

-0.719047 

-0.732727 

-1.028913 

0 

2 

0.183191 

-0.291969 

-0.275917 

-0.362931 

0 

3 

0.081475 

-0.117703 

-0.099753 

-0.118659 

l 

0 

0. 

0.090094 

0.161513 

0.304195 

l 

1 

0. 

0.047004 

0.081646 

0.148907 

l 

2 

0. 

0.0250925 

0.0421838 

0.0743690 

2 

0 

0. 

0.0091060 

0.0124670 

0.0163852 

2 

1 

0. 

0.0053384 

0.0071132 

0.0090464 

3 

0 

0. 

0.0018101 

0.0021753 

0.0023946 


Table 3: Hadronic moments without Eg cuts. 

masses = 4.6 GeV and m c = 1.15 GeV. These numerical values do not correspond to the 
actual b and c pole masses, but they simplify the comparison with [18, 31]. 

In table 1 we present the results of L^ 1 for n = 0,1, 2, 3 and j = 0,1, 2, 3 when no lower cut 
is applied on the charged lepton energy. The first row refers to the total width: we reproduce 
the results of [12, 21, 22] for the four entries. The second to fourth rows refer to the first few 
moments of the leptonic energy; we have checked that we reproduce with good accuracy the 
results of [15, 23] for the first two entries, while the third and fourth entries have not been 
reported so far. 

Table 2 contains the same leptonic moments as table 1, but a 1 GeV lower cut on the lepton 
energy is applied. Here the second entry in each row confirms results by [23], while the third 
and fourth entries are new. 

Reference values for the hadronic moments are given in table 3 for the case without lepton 
energy cut, and in table 4 for Eg > 1 GeV. The values of tables 3 and 4 reproduce the tree-level 
and one-loop results presented in refs. [13] and [16], taking into account the different input 
values. The same holds for a few entries of the third column (second-order BLM) given in [16]. 
We always find acceptable numerical agreement with [18]. More difficult is a comparison with 
the one-loop numerical results presented in [17], as they refer to the physical hadronic invariant 
mass moments, and are expressed in the IS b mass scheme. 
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i 

j 

H ^ 

l_l 

Fv 

l_l 

h (2) 

'll 

IJ 

IT 

0 

0.814811 

-1.439726 

-1.600023 

-2.407833 

0 

1 

0.334079 

-0.577255 

-0.613259 

-0.883227 

0 

2 

0.141115 

-0.234529 

-0.235029 

-0.320474 

0 

3 

0.061196 

-0.095945 

-0.089096 

-0.112775 

1 

0 

0 . 

0.057255 

0.108008 

0.212635 

1 

1 

0 . 

0.028504 

0.052497 

0.100739 

1 

2 

0 . 

0.0144748 

0.0260422 

0.0486621 

2 

0 

0 . 

0.0044430 

0.0066138 

0.0093668 

2 

1 

0 . 

0.0024224 

0.0035595 

0.0048962 

3 

0 

0 . 

0.0006414 

0.0008931 

0.0010505 


Table 4: Hadronic moments with Eg > 1 GeV. 


5 Summary 

We have performed a calculation of the 0(a s ) and 0(o^/3q~ 1 ) corrections to the differential 
rate for inclusive semileptonic b decays. We have calculated all five structure functions; only 
three of them contribute to the decays into massless leptons, while they are all relevant for 
B —> X c tu. Our one-loop calculation checks several existing results, while the BLM calculation 
is far more general than those present in the literature. Since we explicitly subtract the soft 
singularity from the gluon emission contribution, our results allow for an accurate numerical 
implementation, which we provide in a FORTRAN code, available upon request. The code 
computes 0 (a s ) and O(cd?/3o _1 ) corrections to the triple differential distribution and integrates 
over the phase space with arbitrary experimental cuts. A few numerical results are presented 
in Section 4 for the case of massless leptons. 

For what concerns the analyses of semileptonic moments, relevant to the determination of 
V cb |, the present work permits the evaluation of the perturbative contributions to hadronic 
moments at the same level as for the leptonic moments, since all 0 (a s ) and 0{o^(5q^ 1 ) contri¬ 
butions can be now computed. Practical applications have been touched upon in [18] and will 
be expanded in a future publication. 
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A Kinematics 

Let us consider the phase space element for the processes 

Kp) +P(Pv) + AT C , 


(A.l) 


where the states X c contain a charm quark with momentum p', plus a collection of n partons 
with momenta ki,..., k n and masses mi,..., m n . It can be decomposed as 

dq 2 dr 2 

d<f> 3 + n {p'lPiiPviP' ,ki) = — —d(f)2(p;q,r)d(j)2(q]pe,PD)d(f) n+ 1 (r-,p',k i ) i (A.2) 

where 


2 ^ 2 ^ 2 

q- < q < q +; 

2/2^2 
r_ < r < r , ; 


q- = {me + m v ) 2 = 0; q 2 + = (m b - m c - rn x 


m r 


(A.3) 


r_ = (m c + mi + ... + m n ) ; r + = (m b - Jq 2 ) . (A.4) 


In the rest frame of the decaying b quark, the direction of the momentum q of the lepton pair 
is irrelevant, and its modulus is fixed for fixed g 2 , r 2 . Therefore, the first two-body phase space 
in eq. (A.2) gives simply an overall factor: 


dMp '«' r) = (2$2i 0 (2$2^ (2jr)Vi){p - 9 “ r) = lm b ’ 


where 


|^2 = A (q 2 ,r 2 ,ml) 


4 ml 


\{x, y : z) = x 2 + y 2 + z 2 — 2 xy — 2 xz — 2 yz. 


(A.5) 


(A.6) 


We now turn to the phase space factor for the lepton pair. For a given value of q 2 , the 
kinematics is completely fixed by the energy of one of the two leptons, for example the charged 
lepton energy Ee\ 


= ( 4 % < 2 ") 4|i(4) < 9 ~ ~^ 

1 Ee dEe d cos 9 r/ 0 


8vr \q-pe\ 


5{q - E e - \q-p e \). 


(A.7) 


where 6 is the angle formed by the directions of q and pe- The delta function can be used to 
perform the angular integration. We obtain 


dh{q^Pe,Pp) = 


1 dEe 
8tt | q\ ' 


with 


cos 6 = 


2 q°E £ - q 2 
2\q\E e 


q = yq 2 + \q\ = 


2 m 2 + q- — r 2 


2 m b 


From the condition | cos 6* | < 1 we get the limits on Ep. 


Ef, < Ee ^ — 


q° ± \q\ 


(A.8) 


(A.9) 


(A.10) 
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Collecting eqs. (A.2,A.5,A.8) we obtain 


d<j)3+n(p;Pe,Pv,p',k) = 


1287r 4 mb 


dE t dq 2 dr 2 do, l+] ir: k n ). 


(A.11) 


For fixed q 2 and Eg, the variable r 2 is in one-to-one correspondence with the neutrino energy 
E p ; 

r 2 = (p - q) 2 = ml + q 2 - 2 m b {Eg + Ep). (A.12) 

Similarly, the invariant mass of the lepton pair q 2 can be traded for the energy Ex c of the 
hadronic system in the final state; they are related through 


q — ip- r) = ml + r - 2 m b E Xc - 
The differential decay width for the processes in eq. (A.l) is given by 


(A.13) 


d 3 T _ 1 1 1 

dq 2 dr 2 dEg 2 2 m b 128n 4 m b 


Y j d(j) n+l {r-,p l ,k 1 ,...,k n ) \M n {p;pe,pp,p',ki)\ 2 . (A.14) 


At leading order in the weak interactions, we have 

X X 

M n (p;pe,pp,p' } ki) = —^— —M^{p] q,p'.k^M^pt.Pv) ~ —M^{p] q,p'^JM^ p^Pp), 


q' 2 - m 2 w 


where Ad^ip; q,p ', kj ) and Ad ^{q; pg, p p ) are the amplitudes for the processes 

b(p)^W*(q) + X c 
W*(q)^£(p e )+9(pp) 


(A-15) 

(A.16) 
(A.17) 


respectively. We define a leptonic tensor U‘ ,v and a hadronic tensor W^ through 

g 2 E w (pe, p D ) = Ai fM (q-,pe,pp)Ai*"(q-,pg,pp) (A.18) 

g' 2 \V cb \ 2 W^(p,q) = Yf d <l>n+i{r;p',k 1 ,...,k n )Ad l l(p;q,p , ,k i )Ai* n l '{p;q,p , ,k i ).(A.19) 


In terms of L^ and W^, the decay rate is given by 

d 3 r 111 g 4 \V cb \ 2 


dq 2 dr 2 dEg 2 2m b l28n A m b m^ v 


Lgw(pe,Pp) W^ip, q) 


(A.20) 


Finally, we compute the phase space for the hadron system in the cases n — 0,1. For n — 0, 
which is the case at order a 3 and for the virtual contribution at order a s , we have simply 


dMnp') = 7x 7 (27r) 4 A (4) (-r - p') = 2vr 5(r 2 - m 2 ). 

{Z7T) Zp o 


(A.21) 
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For n — 1 we have a charm quark with momentum p' and a gluon with momentum k in the final 
state. Choosing the rest frame of the charm-gluon system, with the 3-axis along the direction 
of p, we find 


k) = (2^ i2n)isW{r ~ p ' ~ k) = pyV 059 - 


(A.22) 


2 A(r 2 , ml, X 2 ) / -> , r 2 — m 2 + A 2 


\k\ 2 = 


ko — y\k\ 2 + A 2 = 


2 V^ 


It will be convenient to introduce the variable 


(A.23) 


t = (p — k) — ml = A — 2(p 0 k 0 — \p\\k\ cosd); dcosd = 


mw 


where 


bf = 


A (q 2 ,r 2 ,m 2 b )_ 0 _ / 2 _ m'i - r 2 ^ <f 


P = Vb1 2 + m fe = 


2 y/r 2 


because of the energy conservation constraint 


Hence 


| p\ 2 + ml = vr 2 + a/ |p| 2 + q 2 . 


t+<t<t-- t± = A 2 - 2p 0 k 0 ± 2\p\\k\. 


(A.24) 


(A.25) 


(A.26) 


(A.27) 
(A.28) 


B Virtual corrections 

In this Appendix, we present in full detail the computation of the order-aq contribution to the 
differential rate due to one-loop corrections to the process in eq. (2.5). We have 


W m^ q) = / (27T) d £ 2 - A 2 \{i + p) 2 - m 2 } [(£ + p 1 ) 2 - m 2 ] 

xTr [(/+ m c )'f p (/ + /+ m c ) 7 M (l —7 5 )(/ + 1> + + m b )Y(l ~ 75 ) 

+ tf> + m b )7 P {f. + i> + m b )7 V {1 - 75) (/ + /+ m c ) 7 P (/+ m c ) 7^(1 - 75) 

-* (Z b + Z c ) W$(p,q), (B.l) 


where 


dE(p,mi) 

dfi /p=mi 


S(/J, rrii) = C F g 


2 r d d i 1 7p(/ + ^ + '^)7 p 

s J (2ir) d £ 2 — A 2 (£ + p) 2 — m 2 ’ 


(B.2) 
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and p' — p — q. With standard techniques, we obtain 

w wv(p,q) = - iC ' F9 's f o dx j Q dy 


(27r) d (£ 2 - A/ 2 ) 3 + . 
2 - d I 

J (2vr) 4 (£ 2 - M 2 ) 3 _ 
’ d 4 £ 8m 2 x(l + x) 

(2<K) d (£ 2 - M 2 ) 2 J 

d d i 2 - d j 

(2tt) 4 (£ 2 - M 2 ) 3 
’ d 4 £ 8m 2 x(l + x) 

(2ir) d (£ 2 - M 2 ) 2 J 

(2vr) 4 (£ 2 - M 2 ) 3 


where we have collected in N\ V all terms quadratic in the loop momentum £, while is 
^-independent. Furthermore, 

M 2 = mix 2 + m 2 c y 2 — uxy + A 2 (l — x — y) (B.3) 

M 2 = m 2 x 2 + A 2 (l — x), (B.4) 

where uj — q 2 — m 2 — rn 2 . Ultraviolet divergences are regulated by dimensional regularization, 
with d — 4 — 2e. Using the properties of symmetric integration, we can replace 

Nr 8 e 2 ^—-^1 L + 4 p^p u + 2ie^ a0 p a q 0 - 2 (q^p u + q v p^ L )} , (B.5) 

d i J 


7vr 8 e v ; \u g^ + 4p^ 

CL 

and hence, by comparison with eq. (2.6), 

N^nSiu) ^8£ : 


,(2 -df 


Nr IT S(u) 8 Wft (p, q)- (B.6) 

Ultraviolet divergences are immediately seen to cancel, and setting d — 4 we are left with 

r< n { t t r 1 r 1 ~ x N^ i 

W^ v (p, q) = -^|(1 + I 0 - 2 J) Wfi(p, q) + - S(u) ^ dx ^ dy - j, (B.7) 


where we have defined 

h = L dx L dyXog MfMj 

After some algebra we find 

AT = -16a; [a; + 4p> v + 2 ie^ a0 p a q 0 - 2{p^q v + j+<f) 

+ 16 luMq + (to 2 — Luml — 2mlm 2 c )x + ( u 2 — a vnr? c — 2ml r m 2 c )y g^ v 

—64 (a; + m 2 c + m\)xy — u(x + y) p^p u 

+32 Ml + (a; - m 2 b )x + (a; - m 2 c )y ie^ a0 p a q 0 

+64m|x(l — y) q^ L q u 

+32 (a; + 2 ml)xy — (u + — (a; — m 2 c )y (p^q u + p v q M ), (B.9) 
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where 


(B.10) 


7i /r2 2 2 i 2 2 

M 0 = m b x + m c y — coxy. 

Our final result for the virtual contribution is therefore 


Wfi y (p,q) = - C F a s {V 0 (q 2 ,X ) [u + AvW + 2ie^v a qy - 2{q‘*v v + gY) 

- lh(<f, A) <T + 0 2 (<f, A) u V + ^(g 2 , A) u a g> 

+ 14(<f, A) <f <f + 0 5 (g 2 , A) (iPq u + g"/)} 5(u), 


and 


(B-ll) 


where 


Vote 2 . A) = i 

(1 + Jo — 2J — 2/i) 



to 

1 

1 

2 

a; A' + (a; 2 — cum 2 — 2 m 2 m; 

l)i 

x + (^ 2 — 


m 

'6 




V 2 {q\ A) = - 

4 

(cn + m 2 c + ml)I xy — u(I x + I y ) 


V 3 (q 2 , A) = 2 

A' + (ca - ml)I x + (ca - m 2 )4 



V 4 (q 2 , A) = 4m 2 (/ 

a? Ixy) 



O s (g 2 ,A) = 2 

(cu + 2 ml)I xy - (u + ml)I x - ( 

[co 

- ™c)4] 


.2 2 \ 


(B.12) 


K = 
h = 
4 = 
4 = 

Ixy = 




/>r^§ 

/•I /*1—# 

L dx l dv M? 

[1 rl-x x 

L dx i dv w 

y 


r 1 


da; 


da; 


r»l—# 


dy 


Ad 2 


r»l—# 


a;a/ 

dy w 


(B.13) 


The integrals Ji and J are divergent for A = 0; this is precisely the soft singularity which is 
needed in order to cancel the analogous divergent terms in the real emission contribution. Also 
in this case, it is convenient to isolate the soft logarithm. To this purpose, we perform the 
change in the integration variables 


x = XiX 2 (B.14) 

y = x 1 (l-x 2 ) (B.15) 

dx dy = X\ dxidx 2 ; 0 < Xi < 1; 0 < x 2 < 1. (B.16) 

In terms of the new variables, we have 

M 2 = P(x 2 )x\ + A 2 (l-a; 1 ), (B.1T) 
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where 



Thus, 



h = 

oj 

f 1 dx 2 f 1 

lo Jo 


OJ 

f 1 dx 2 


2 

Jo P(x 2 ) 


c 0 

r 1 dx 2 


2 

•JO P(x 2 ) 


P(x 2 ) = mix 2 + m 2 (l - £ 2 ) 2 - ^ 2(1 - z 2 ). 


(B.18) 


xi 


0 ' P(x 2 ) x\ + A 2 (l - xi) 


dx 1 


log- 


A 2 


2P(x 2 )xi — A 2 


+ 


+ 


cixi 


a; f 1 dx 2 f 1 


A 2 


cixi 


2 Jo P(x 2 ) Jo ' P(x 2 ) x 2 + A 2 (l — xi) ’ 


(B.19) 


The coefficient of log A 2 can be computed explicitly. We find 


u> r 1 dx 2 1 1 + a 

2 Jo P(x 2 ) 2a ® 1 — a 


(B.20) 


where a = Hence, 


/, = — log log — 

2 a 1 — a m\,m c 

u r 1 dx 2 , P(x 2 ) ca /•> dx 2 r 1 , 

log- + - / T77 T / <fcl 


A 2 


2 Jo P(x 2 ) rnb r m c 2 Jo P(x 2 ) Jo P(x 2 ) xf + A 2 (l — aq) 


(B.21) 


Similarly, 


nl m 2 a;(l — a; 2 ) 


1 /•' , 2m?x — A 2 
1 dx ■ 


+ 


1 r 1 , A 2 — 2mix 3 
1 dx ■ 


1 0 m 2 a; 2 + A 2 (l — x) 2 Jo m 2 a; 2 + A 2 (l — x) 2 Jo m 2 a; 2 + A 2 (l — x) 

1 m? + i[' dx , (B.22) 


2 l0g A 2 ' 2 Jo "~m|;x 2 + A 2 (l -z) 


and therefore 


J = - log —^-f J f 1 dx 

nii/nic 2 Jo 


A 2 — 2 mix 3 


+ 


A 2 — 2mJx 3 


mix 2 + A 2 (l — x) m 2 a; 2 + A 2 (l — x) 


(B.23) 


C Analytic expressions for one-loop moments 

We present here a few examples of 0(a s ) corrections to the moments of u 1 Eq at arbitrary values 
of the lower cut on the charged lepton energy. Keeping in mind eqs. (3.2) and (4.2) we define 

]\/r r 0 1 1 ) 

^ F u 


with T u = G 2 F \V c b\ 2 ml/(192Tc 3 ), £ = 2E‘j, nt /m^ and L^ = ln(l — £). 
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The moments with % > 0 receive contributions from bremsstrahlung diagrams only and are 
significantly easier to express in analytic form. Here are the lowest moments: 


_ —2873p 5 , p 4 (l-62C +108e 2 -93e 3 ) , (1 - O' (91 - 298e - 52£ 2 + 90) 

M 10 — 7WXX-1-V - T“-V)-T 


1800 


24(1 - a 


600 


p 3 (i7 +109 e - 83 e 2 + is a p 2 (-79 - 438 e +222 e 2 + 25 e 3 ) 


6-ee 

p(i-o (7i + io7e+45e 2 -5ia 

24 


18 


+ 


p 3 (63 - 37e~14a 
3 - 3 ^ 


p 2 (-146 + 195 e - 60 e 2 + 2 a (1 - O 2 (96 + 17£ - 62 £ 2 + 9 a 


6 


120 


p 4 (-60 + me - use 2 + 17 a + p (9-2 oe+iso - 7 a 


+ 


24 (1 — e) z 

pH 


4p 3 (3-o-^+p 2 (s-oe+e 3 ) 


L^ In 


6 

i-e 

p 


+ 


7p 5 p 3 (i + ne-ioa 


15 


p 4 5-ioe + 6e 2 2 5 r , 2 „ . 

—- 9 7 + p - - 8e+ 5 e - — + 

2(1-0 V 3 6 / 


e 3 \ , 3p ( 1 - 20 +a 


i-e 

lnp 


(C.l) 


6341 p 6 , p 2 (l-0 (-101-20ie + 6ie 2 + 4ia , p 3 (43 + 528e - 2850 + 47a 

“ 10800“ + 48 + 54 
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